Spectrum and Wave Functions of Excited States in Lattice Gauge Theory 
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We suggest a new method to compute the spectrum and wave functions of excited states. We 
construct a stochastic basis of Bargmann link states, drawn from a physical probability density 
distribution and compute transition amplitudes between stochastic basis states. From such transi- 
tion matrix we extract wave functions and the energy spectrum. We apply this method to U (1)2+1 
lattice gauge theory. As a test we compute the energy spectrum, wave functions and thermodynam- 
ical functions of the electric Hamiltonian and compare it with analytical results. We find excellent 
agreement. We observe scaling of energies and wave functions in the variable of time. We also 
present first results on a small lattice for the full Hamiltonian including the magnetic term. 



INTRODUCTION 



Much progress has been made in lattice gauge the- 
ory on the computation of hadron spectra jl|, H, 0, 
0, M 0| ■ Such spectra can be extracted from a matrix 
of 2-point correlation functions built from extended 
operators with suitable quantum numbers applied to 
the ground state [3] . Here we present a new approach 
for the computation of spectra and wave functions 
suggesting to compute transition matrix elements be- 
tween states chosen from a stochastic basis with re- 
spect to gauge degrees of freedom. We consider tran- 
sition amplitudes 

M li , v (T) = (r„\exp[-HT/h]\? li ) , fi,u = 1,...,N , 

(1) 

where \T V ) denotes a time-independent Bargmann- 
link state, i.e., a configuration of link variables f7y 
assigned to all of the links ij on the spatial lattice. It 
is crucial to choose states |T„) which are physically 
relevant and important. We use as stochastic tech- 
nique Monte Carlo with importance sampling to sam- 
ple states from a large variety of possibilities. The 
stochastic basis states are closely related to equilib- 
rium path configurations in the Euclidean path inte- 
gral. From such matrix M^ V (T) we extract a spec- 
trum and wave functions of an effective Hamiltonian 
- the so called Monte Carlo Hamiltonian - being valid 
in a low energy, respectively low-temperature window. 
The Monte Carlo Hamiltonian has been suggested 
1999 [§]. In field theory, the Monte Carlo Hamilto- 
nian has been applied to the 1 + 1 Klein-Gordon model 
for the computation of the spectrum and thermody- 
namical functions 1(J II , 12|, HI 1 an d likewise to the 



1 + 1 scalar model for the computation of th e sp ec- 
trum and thermodynamical functions [l4, IB], Hf- A 



first step towards the Monte Carlo Hamiltonian in lat- 
tice gauge theory has been made in (l7l | by computing 
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transition amplitudes of f7(l) gauge theory. Here we 
construct the Monte Carlo Hamiltonian for {7(1)2+1 
lattice gauge theory and apply it to compute the en- 
ergy spectrum of excited states, the corresponding 
wave functions and thermodynamical functions. Here 
we want to show the working of the method applied 
to an Abelian, but non-trivial model in lattice gauge 
theory. In the case of the electric Hamiltonian we con- 
sider spatial lattice volumes up to 10 2 , while for the 
full Hamiltonan we present only preliminary results 
on a 2 2 lattice. Because of the smallness of the lat- 
tices we do not consider here the quantum continuum 
limit (o — > 0). Rather, we try to give a careful anal- 
ysis of the origin and size of errors, which determine 
the limitations of the method. 



{7(l) 2+ i LATTICE GAUGE THEORY 



The lattice action is given by [18 1 
1 a 



S[U] 



+ E [l-^n)] 

U space-like 

= S e l ec [U] + S ma g n [U] . 

The corresponding lattice Hamiltonian is given by [I 

- 2 — - 1 



(2) 



H 



■L 

2a 



g 2 a 



£ [l-Re(U n )] 



(3) 



representing the electric term and a magnetic term, 
respectively. The electric term is built from the oper- 
ator lij which represents the electric flux strings. Its 
eigenstates are 



hj\Xij >= Ay I > , Ay = 0, ±1, ±2, . . . . (4) 
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FIG. 1: Scaling of energy levels from the Monte Carlo 
Hamiltonian. Comparison of electric Hamiltonian (left) 
with full Hamiltonian (right). The bars represent statisti- 
cal errors. 2 2 spatial lattice, a s = at = 1, g = 1, N = 80. 



For each link ij, the states |A) form a complete or- 
thogonal basis, 

J2 |A><A| = 1 , (X'\X) = 8 KX . (5) 

A=0,±l,±2,... 

The magnetic term is built from link operators Uij . It 
has the eigenstates 



U ij \U^) = U ij p^) 



(6) 



For each link ij, this basis is a complete orthogonal 
basis, 



dU \U)(U\ = 1 , (U'\U) = 5{U' - U) . (7) 



The transition from the flux string basis to the link 
basis is determined by the commutator [I, U] = — U 
and yields [17 1 



We construct Euclidean transition matrix elements 
M Uf>Ui = (U f \iLexp[-HT\\Ui) 

U f ,T 

[dU] exp[-S[U\] 

Ui,0 

(9) 

where the operator II (commuting with the Hamilto- 
nian) denotes a projection operator states onto gauge 
invariant states. 
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FIG. 2: Electric part of Hamiltonian. Comparison of 
exact results with those from Monte Carlo Hamiltonian. 
Left: Low-lying spectrum. 10 2 spatial lattice, a s = a t — 1, 
g — 1.5, T — 4.2, N = 1000. Right: Specific heat C(/3), 2 2 
spatial lattice, a s = at = 1, g = 1, N = 500. Error bars 
indicate statistical errors. 



STOCHASTIC BASIS 

For a system with a given Hamiltonian H the phys- 
ically motivated choice for the distribution is given by 
the transition amplitude in imaginary time, involving 
the same Hamiltonian, 



< X\U >= (U) x . 



(8) 



P{U) = (U\exp[-HT/h}\U imt ) , 



(10) 
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where Uinit is some suitably chosen fixed spatial lat- 
tice configuration (Bargmann state). This function 
is suitable as probability distribution because it is a 
positive function P(U) > 0. Physically relevant con- 
figurations (Bargmann states) can be drawn from this 
distribution by expression P(U) as path integral and 
doing the sampling via Monte Carlo. As alternative 
one may consider the distribution given by the tran- 
sition amplitude from the electric Hamiltonian, 

P{U) = (U\exv[-H elec T/h}\U init ) , (11) 

which is an analytically computable function. Both 
of the above distributions involve a time parameter 
T, which determines the "width" of the distribution. 
Such time parameter needs to be tuned. As a general 
rule, we used to choose T such that the energy spec- 
trum falls into a scaling window of eigenvalues (see 
below) . 



TRANSITION AMPLITUDES FROM 
ELECTRIC HAMILTONIAN 

Applying the Peter- Weyl theorem [2(| to the group 
U(l) allows to expand link states in terms of irre- 
ducible representation matrices. In the case of the 
group U(l), the Peter- Weyl theorem is equivalent to 
Fourier expansion [2(| . The Peter- Weyl theorem holds 
more generally for groups SU(N) [21(. Applying the 
Peter- Weyl theorem for the group U(l) to the transi- 
tion amplitude between two single link states reads 



< U fi \exp[-H e i ec T/h]\U in > 
g*hT 



E 



exp[ 



n=0,±l,±2, 



2a 



n ] cos[n(aj„ - a/,)] 



(12) 



In mathematical terms n (running over 0, ±1, ±2, . . . ) 
denotes the index of the irreducible representation. 
(U) n denotes the irreducible representation of group 
element U with representation index (quantum num- 
ber) n. In physical terms, n represents the number 
of electric flux lines. The Hamiltonian H e i ec is a 
Casimir, which is diagonal in the representation in- 
dex n. The link variables have been parametrized via 
U = cxp[ia]. 

The Peter- Weyl theorem is also useful for the con- 
struction of gauge invariant states. For example, let 
us take a spatial lattice consisting of four links or- 
dered to form a plaquette and consider the transition 
amplitude between initial and final Bargmann states. 
In order to make the amplitude gauge invariant, we 
carry out the group integral over the gauge orbit at 
each node. The group integral generates Gauss' law 
enforcing conservation of the number of flux lines at 
each vertex. By defining the plaquette angle 9 p i aq = 
"12 + "23 + "34 + "4i and A6 p i aq = 9£ aq ~ 6 l p n laq , we 



obtain the final expression of the gauge invariant am- 
plitude, 



(U fi \ileM-HeiecT/h)\U 
9 2 hT 



J2 ex P 



ra=0,±l,±2,. 



2a 



4n< 



cos 



nA9 p i aQ 



(13) 



Here n denotes the number of closed plaquette loops. 
The result is built from plaquettes which are closed 
loops of consecutive link variables forming the smallest 
non-local gauge invariant objects on the lattice. The 
eigenvalue of the electric field E 2 corresponds to the 
contribution from n plaquette loops (on top of each 
other). The result only depends on the number of 
plaquette loops and the difference between initial and 
final plaquette angles. 



TEST OF MONTE CARLO HAMILTONIAN 

A. Electric Hamiltonian. In order to test the Monte 
Carlo Hamiltonian, we first consider the electric part 
of the Hamiltonian. This is a good test bed, because 
the spectrum of the electric Hamiltonian can be com- 
puted analytically. Fig. pQ (left) shows the low- lying 
part of the spectrum as function of the time param- 
eter (3 = T occuring in the matrix elements. The 
physical spectrum should be independent of the time 
parameter. In the numerical results this is reflected by 
the existence of scaling windows (region of flat line). 
The size of such scaling window depends on the par- 
ticular energy-level and decreases with increasing en- 
ergy. We found that the size of the scaling window S n 
can be described approximately by an exponential law 
S n oc exp[—crE n ]. Such behavior of decreasing scaling 
windows can be understood from the property that 
exp[— HT/H] projects onto the ground state for large 
T (Feynman-Kac theorem). Higher levels become ex- 
ponentially suppressed by the dominant ground state 
and can survive only for short times T. We looked also 
for scaling windows in the corresponding wave func- 
tions. In particular, we have studied (e^|$„), i.e. the 
expansion coefficient of wave function $„ in terms of 
the stochastic basis function e M . Such expansion co- 
efficients for the first energy levels expanded in terms 
of the first basis function also display scaling windows 
(not shown). Like the size of the energy scaling win- 
dow decreases with increasing energy E n , also the size 
of the wave function scaling window decreases with the 
level index n of energy. 

For the case of low-lying spectrum using a lattice 
of size 10 2 a comparison of results from the Monte 
Carlo Hamiltonian with the exact spectrum is shown 
in Fig. [2] (left). The high degeneracy is due to 1- and 
2-plaquctte states located anywhere on the 10 2 lattice 
(such degeneracy will be lifted when taking the mag- 
netic term into account). The figure shows, firstly, 
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that the Monte Carlo Hamiltonian captures almost 
all of the degenerate states and secondly, reproduces 
the exact energies with small error. The quality of 
the energy spectrum of the Monte Carlo Hamiltonian 
can be seen also from a look at thermodynamical 
functions. In the case of the electric Hamiltonian, 
the energy spectrum Eq, E± , . . . can be computed 
analytically. Thermodynamical functions can be 
expressed via those energies. We have computed 
average energy, free energy, entropy and specific 
heat, and compared the results from the Monte Carlo 
Hamiltonian with the exact one. The result for the 
specific heat is displayed in Fig. [2] (right). In general, 
one observes very good agreement in the regime 
of large f3, i.e., the low temperature regime. For 
small values of (3 some disagreement becomes visible, 
reflecting the fact that the precision of higher enegy 
levels of the Monte Carlo Hamiltonian is limited 
(their scaling windows go to zero). 

B. Including the magnetic term: full Hamiltonian. 
Finally, we consider the gauge invariant transition am- 
plitude under the full Hamiltonian (Eq[T]). Although 
it can be expressed in terms of a path integral with 
the lattice action (Wilson action), this is numerically 
not suitable, because Monte Carlo with importance 
sampling only allows to compute ratios of transition 
amplitudes. Hence, we factorize the above amplitude 
into two terms, one being analytically computable and 
the other one being given by the ratio of transition 
amplitudes computable via Monte Carlo, 

M^(T) = (Cgn exp[-H elec T/h]\U v ) 

x (tgn cxp[-gr//t]|^) 

(ign exp[~H elec T/h]\U u ) ' 

Taking into account the magnetic term allows to ob- 
tain the full Monte Carlo Hamiltonian. First results 
on scaling of its low-lying energy spectrum are shown 
in Fig. [1] (right). These results correspond to a small 
lattice (2 2 ) and also a small number of basis func- 
tions. The results show scaling windows for the low- 
est five energy levels. Compared to the scaling behav- 
ior observed in the electric Hamiltonian (Fig. [T] left) 
fewer levels show scaling, and the scaling windows are 
smaller. This can be understood from the fact that 
the ratio of matrix elements in Eq. (| 14|) has been de- 
termined via Monte Carlo from path integrals, which 
carries statistical errors in the order of a few percent. 
From this observation we conclude that the numerical 
resolution of energy levels and the size of scaling win- 
dows of the full Hamiltonian is essentially determined 
by the statistical error occuring in the numerical cal- 
culation of the ratio of matrix elements. Results with 
better statistics, a larger stochastic basis and larger 
lattice volumes are under way. 
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